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In this work, we have considered Cardassian Universe in Hoi'ava-Lifshitz gravity. Four types of 
Cardassian Universe models i.e., polytropic/power law, modified polytropic, exponential and mod- 
. ified exponential models have been considered for accelerating models. The natures of statefinder 

parameters, deceleration parameter, Om diagnostic and EoS parameters have been investigated for 
all types of Cardassian models in Hofava-Lifshitz gravity. 

^ ■ I. INTRODUCTION 



Recent observations of the luminosity of type la supernovae indicate an accelerated expansion of the universe 
' [1-7] and lead to the search for a new type of matter, called dark energy, which violates the strong energy condi- 
tion [8-10]. For explaining the accelerating expansion of the universe, large number of cosmological models have 
f~| ■ been proposed. Dark energy model is proposed by assuming an energy component with negative pressure in the 
D ' universe, this dark energy dominates the total energy density of the universe and drives its acceleration of ex- 
pansion at recent times. Some authors have explored possible explanations for the acceleration: (i) cosmological 
c/J ' constant, (ii) quintessence [11-17], (iii) gravitational leakage into extra dimensions [18-19]. Accelerated expan- 
^ , sion also depends on other observations such as the cosmic microwave background (CMB) [20] and galaxy power 
c/3 ■ spectra [21]. Another reason of acceleration is the geometric effect of the general relativity fails in the present 
cosmic time space scale. The model proposed by Freese and Lewisan, called Cardassian model [22-26], assume 
that the universe is flat and accelerating which consists only of matter and radiation [27]. The geometry is flat 
as required by measurements of the cosmic background radiation, so that there are no curvature terms in the 
equation and no vacuum term in the equation and so the model does not address the cosmological constant. In 
Cardassian Model, we take g{p) as a function of p that returns simply to Pp/3 at early epochs, where P = SnGc. 



On the basis of the recent observation one can state that if Einstein's theory of gravity is acceptable on 
cosmological scales, then our universe must dominated by a mysterious form of energy called dark energy. 
Recently Hofava [28-29] proposed a new theory of gravity, which renormalizable with higher spatial derivatives 
in four dimensions which reduces to Einstein's gravity with non vanishing cosmological constant in IR but 
with improved UV behaviours. In Lifshitz [30] scalar field theory the time dimension has weight 3 if a space 
dimension has weight 1 and this theory is called Hofava-Lifshitz gravitational theory. Hofava-Lifshitz gravity 
has been studied and extended in detail and applied as a cosmological framework of the universe [31-33]. 

. , — I , As so many cosmological models have been developed, so for discrimination between these contenders, Sahni 
^ ■ et al [34] proposed a new geometrical diagnostic named the statefinder pair {r, s}, where r is generated from 
?H [ the scalar factor a and its derivatives with respect to the cosmic time t and s is a simple combination of r 
. . . 1 and the deceleration parameter q. Clear differences for the evolutionary trajectories in the r-s plane have been 
found. In this work, we have discussed four different types of Cardassian Universe models in Hofava-Lifshitz 
gravity. In every case we find the statefinder parameters, deceleration parameter and Om diagnostic [34-36]. 



II. CARDASSIAN UNIVERSE IN HORAVA-LIFSHITZ GRAVITY 

The Arnowitt-Deser-Misner formalism of the full metric is written as [37], 
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ds^ = -N^df + g^j{dx' + N'dt){dx^ + NHt) (1) 
Under the detailed balance condition the full action condition of Hoi^ava-Lifshitz gravity is given by, 



S = I dtd^x^N 



8(3A- 1) 



1 -4A 



-R^ + AR- 3A^ 



(2) 



where the extrinsic curvature and Cotton tensor is defined as, Kij ~ jj^idij ~ ^i^j ^ ^ j^i) o,nd C"-' = 
^^Vfc(i?^ — jRSj). The covariant derivatives are defined w.r.t. the spatial metric gij. e^^^ is the totally 
antisymmetric unit tensor, A is a dimensionless coupling constant and the variable ft , uj and /i are constants 
with mass dimensions — f, 0, f respectively. Also A is a positive constant, which as usual is related to the 
cosmological constant in the IR limit. 

Now, in order to focus on cosmological frameworks, we impose the so called projectability condition and 
use a Friedmann-Robertson- Walker (FRW) metric we get, N = l^gij = a'^{t)%j,N^ — with ^ijdx^dx^ = 

■y'^^^-i + r^d^i^ where fc = 0, —1, +1 corresponding to flat, open and closed respectively. By varying N and g^, 
we obtain the non- vanishing equations of motions: 



and 



H +-W 



6(3A-f) 6(3A-f) 



3KVfc' 3k2^2a2 



8(3A-f)a4 8(3A-f) 



4(3A-f)'^ 4(3A-f) 



3kVA^ 

8(3A- f)a4 ~ 8(3A- f) 



8(3A- f)2a2 



f6(3A-l)2a2 



(3) 



(4) 



where H = ^ \s the Hubble parameter. Here Gc and G is defined as, Gc = 16^(3^-1) ^^"^ ^ = Mt: where Gc 
is the "cosmological" Newton's constant and G is the "gravitational" Newton's constant. 



We can re- write the above equations as [38], 



and 



2- _fc^ A 

3 ' 2 iW* ^ 2 



Pp fc2 3A 
~2~ " 4A^ ^ T 



(5) 



(6) 



Here, p and p are respectively the energy density and pressure of the universe, P = 8itGc and choosing 
87rG = 1. 

Freese and Lewis [22] constructed Cardassian universe models so that, in Cardassian models the universe 
is flat and accelerating, and yet contains only matter (baryonic or not) and radiation. The above equation 
governing the expansion of the universe is modified to. 



A 

2Aa4 + 2 



(7) 



which gives. 



H 



k Pq{p) fc2 A 
1 ilLi. A \ 



a 



2Aa4 



(8) 
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where p is the total energy density of matter and radiation and we will neglect the contribution of radiation 
for the late-time evolution of the universe. The function g{p) reduces to p in the early universe. Now, 

g{p) = p,n + Pc and p = p,n + Pc (9) 

where pm and pc are the energy densities of matter and Cardassian term of the universe respectively and pm 
and Pc arc the pressure of matter and Cardassian term of the universe respectively. For any suitable Cardassian 
model, the following requirements that should be fulfilled, (i) The function g(p) should returns to the usual 
form of p at early epochs in order to recover the thermal history of the standard cosmological. (ii) g{p) should 
takes a different form at late times when z ~ 0(1) in order to drive an accelerated expansion as indicated by 
the observation of SNela. (iii) The classical solution of the expansion should be stable, and the sound speed 
of classical perturbations of the total cosmological fluid around homogeneous FRW solutions cannot be negative. 



In order to guarantee the classical solution of the expansion is stable, the sound speed cj. of classical pertur- 
bations of the total cosmological fluid around homogeneous FRW solutions should always be greater than zero. 
If the expansion of the universe is adiabatic, the sound speed of total cosmological fluid can be represented by. 

Now the matter conservation equation gives, 

Prn + 3H{p,n + Pm) = (11) 

and total fluid conservation equation gives, 

g(^) + 3i/(5(p)+p) = (12) 



From the above two equations we get. 



Pc + ^H{pc+Pc)=Q (13) 



From where we get. 



and 



Pc = (Pm + Pni)^i^ ~ dip) - Pm (14) 
Opm 



P = Pm + Pc = {Pm +Pm)^i^ - g{p) = Pm(l + Wm)^^^ ^ g{p) (15) 

apm OPm 

Now for dark matter ~ WmPm, combining this with (11) we get, 

= Poa~'''+"'"^ and p^ = pow^a-^^^+^'-^ (16) 

where po be the integrating constant. Now the relation between scale factor a and the redshift z is given by 
a = 1/(1 -|- z) i.e. z = (1/a) — 1, and we replace all a by z, where value of z is taken z > —1. 



III. STATEFINDER DIAGNOSTICS, DECELERATION PARAMETER AND Om 

DIAGNOSTIC 

The flat Friedmann model which is analyzed in terms of the statefinder parameters. The trajectories in 
the {r, s} plane of different cosmological models shows different behavior. The statefinder diagnostic of SNAP 
observations used to discriminate between different dark energy models. The statefinder diagnostic pair is 
constructed from the scale factor a{t). The statefinder diagnostic pair is denoted as {r,s} and defined as [34], 



(17) 
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where q is the deceleration parameter given by, 9 = ^ff • If universe is present with dark matter then the 
parameters can be expressed as, 

r = l+ I , f |^(Pc +Pc) + (18) 

2(Pc+Pm) \OPc Opm J 

^-^(Pc+Pc) + ^f^{Pm+Prn)] (19) 



{Pc +Prn) \dpc dp, 

and 



'^-UU'-^] (20) 



2 2VPc + A 

As a complementary to {t", s}, a new diagnostic called Om has been recently proposed, which helps to 
distinguish the present matter density contrast in different models more effectively. The new diagnostic of dark 
energy Om is introduced to differentiate ACDM from other DE models. Om diagnostic is defined as [35], 



Omi^) = (21) 

Thus Om involves only the first derivative of the scale factor through the Hubble parameter and is easier to 
reconstruct from observational data. Om is a constant in ACDM model, since it is independent of redshift z and 
it provides a null test of cosmological constant. Om diagnostic can distinguish DE models with less dependence 
on matter density relative to the EOS of DE. 

IV. DIFFERENT MODELS OF CARDASSIAN UNIVERSE 

A. Polytropic/Power Law Model (PL) 

The simplest model is the power law (PL) Cardassian model where g{p) = p + Bp""^, with B and n < 2/3 are 
two constants and the additional term p" satisfies many observational constraints such as if the first Dopplcr 
peak of the CMB is slightly shifted, the universe is rather older, and the early structure formation z > 1 is 
unaffected. In this model g{p) is defined as [23], 



gip) = Pn 



Pcard 



(22) 



where pcard is a characteristic constant energy density and n is a dimensionless positive constants. So from (9), 
(14), (16) and (22) we get, 

_ pg(l + ^)3»(i+.«„) 
Pc - (^TTI) 1^<^J 

Pcard 

and 

_ (n + nu;,„-l)pg(l + z)3"(i+^«^) 

P card 



So the equation of states are given by. 



Wc = — = {n + nwm — 1) (25) 

Pc 




FIG. 4: The variation of s against r for Wm ~ -01, n = .5, Pcard = 



.1, A = .01, po = .001, k = l,l = .9,Ho^ 72. 



6 




FIG. 5: The variation of q against H for Wm ~ .01, n = .5, pcard = .1, A = .01, po = .001, k — 1,1 — .9, Ho ~ 72. 
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FIG. 6: The variation of Om against H for Wm ~ -01, ?i = .5, pcard = .1, A = .01, po = .001, fc = 1, / = .9, Hq — 72. 



and 



Wm + {n + nwm - 1) 



\ pcard 



9{p) 



1 



po(l+^)-^(i+'^"') 

Pcard 



(n-1) 



(26) 



From above, we see that Wc is a constant. From (16), (18), (19), (20), (23) and (24) we get the statefinder 
parameters and deceleration parameter. 



(1 + z)3po(2 + 9u;„,(l + w^)) + (1 + z)-3--" (iilfl! 



'Po 



) Pcard(2 + 9n(l + w;m)(-l + ^^ + nWm)) 



2(l + 2)-Vo + 2(l + z)-3«'" ( 



(l + 2)3(l + ^n,)po 



) 



Pcard 



and 



(1 + z)Vo«^m(ll + 9u;„0 + (1 + z)-3-". 


r(i+ 




n 

Pcard(-1 + n + nw™)(2 + 9n(l + w„)) 


pcard ) 


2(1 + zfp^tUrr^ + 2(1 + 2)-3™™ 


Pcard 




n 

Pcard (-1 +n + nWra) 



(27) 



(28) 



(1 + z) Vo(l + 3u;„,) + (1 + z)-3"''" (ii±£)!!ilI:ilBl) " p,,,<j(-2 + 3n(l + w^)) 
2(1 + z)3po + 2(1 + z)-3»™ ( (l±£)fi!;pi£o 

\ ycard I 



(29) 



Fig.l represents the variation of g against i7. Fig. 2 represents the variation of w against H. From (9), (10), 
(16), (23) and (24) we plot the graph which shows the variation of the square of velocity of sound against H 
is given in Fig. 3. Fig. 4 represents the variation of s against r and Fig. 5 represents the variation of q against 
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H. From (8), (21) and (22) we get the value of Om. Fig. 6 represents the variation of Om against H . The 
values are taken as, Wm = -01, n = .5, Pcard = .1, A = .01, po = .001, fc = 1, Z = .9, Hq = 72. From the figures, we 
see that g,w,q and Om decrease as H decreases and Hcs between and 1. The parameter s increases and 
keeps positive sign as r decreases from positive to negative values. 



B. Modified Polytropic Model (MP) 

Modified polytropic Cardassian model is a sUght modification of the previous one, can be used on all scales, 
but it docs not quite fit the criteria of the Cardassian model as defined above. At late times in the future of 
the Universe, when /9,„ <C Pcard, this model becomes cosmological constant dominated with A = Pcard- This 
energy density is very similar to a model which motivated by gravitational leakage into extra dimensions. In 
this model g{p) is defined as [25], 



gip) = Pr, 



Pcard 



i(n-l)' 



(30) 



where Pcard is a characteristic constant energy density and a 7^ and n are dimensionlcss positive constants. 
So from (9), (14), (16) and (30) we get. 



Pc = il + z)3(i+^"'"Vo(-l + (1 + XMP)i) 



(31) 



and 



Pc^{l + z)3(i+™'")po(XMp(l + Xmp)-'+° - + (1 + Xmp)-'+° {w,n - XMp{n + (-2 + n)w,n))) (32) 
where. 



X 



MP 



(1 + z)-3(i+"'")p.„.d 



Po 



0.0 0.5 1.0 1.5 2.0 2.5 



FIG. 7: The variation of g against H for a = .5, Wm = .01, n = .5, pcard = .1, A = .01, po = .001, k = 1,1 = .9, Ho = 72 
So the equation of state is given by, 

_ (1 + z)3(i+-'")po(^Mp(l + XMp)-^+i - w,n + (1 + Xmp)-^+^{w„, - XMp{n + (-2 + ?i)w„0)) 



(1 + ;^)3(l + u,„Opo(-l + (1 + Xmp)^) 



(33) 



and 



P0Wi 



^(1 + ^)3(l+«.™) + (1 + z)3(l+-™)po(^Afp(l + Xmp)-^ + ^ - W,n + (1 + ^Mp)-^+ « (u;,„ - ^Mp(n + (-2 + n)w„,))) 



po(l + ^)3(1+-''") + (1 + z)3(l+-")po(-l + (1 + XAfp)i) 



(34) 



8 




FIG. 8: The variation of Wc against H for a — .5, Wm = .01, n = .5, pcard = .1, A = .01, po = .001, k — 1,1 = .9, Ho = 72. 




FIG. 9: The variation of w against H for a — .5, Wm = .01, n — .5, pcard = .1, A = .01, pa = .001, k = 1,1 — .9, Ha = 72. 



I.O 
0.9 
0.8 



0.7 
^ 0.6 
0.5 
0.4 



0.0 



0.5 



1.0 1.5 
H 



2.0 



2.5 



FIG. 10: The variation of against H for a — .5, Wm ~ .01, n = .5, pcard = .1, A = .01, pa = .001, k — 1,1 = .9, Ho = 72. 
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FIG. 11: The variation of s against r for a — .5, Wm = .01, n — .5, pcard = .1, A = .01, po = .001, k = 1,1 = .9, Ho = 72. 



From (16), (18), (19), (20), (31) and (32) we get the statefinder parameters and deceleration parameter, 
r= [2 + 9w,n + 9wl^ + XMpi-5 - 6iy,„ + 3?i(l + w,„))(-4 - 6w„^ + 3n(l + uv„)) + Xmp(13 + 9a(l + w,nf 



H 



FIG. 12: The variation of q against H for a — .5, Wm ~ .01, n = .5, pcard = .1, A = .01, po = .001, k = 1,1 = .9, Ho = 72. 



FIG. 13: The variation of Om against H for a = .5, Wm ~ .01, n = .5, pcard = .1, A = .01, po = .001, k — 1,1 = .9, Ho = 72. 



+971^0(1 + w,„)^ + 9w,„(5 + 4w,„) - 9n(l + w„0(l + 2w„i + 2a(l + w,„)))] / [2(1 + XmpY 



(35) 



and 



s = — 
2 



20 — 9n + 18i(;„i — 9nwm — 



9(-1 + ?i)(-1 + q)(1 + w„0 
1 + Xmp 



9(-1 + n)awm{l + Wm) 
-Wm + Xmp{-1 + n + (-2 + n)wm) 



(36) 



^=2 



4 — Sn + 6k7„ 



3m«„ 



3(-l + n)(l + w„,) 
1 + Xmp 



(37) 



Fig. 7, 8 and 9 represents the variation of g, Wc and w against H respectively. From (9), (10), (16), (31) 
and (32) we plot the graph which shows the variation of the square of velocity of sound against H is given 
in Fig. 10, Fig. 11 represents the variation of s against r and Fig. 12 represents the variation of q against H. 
From (8), (21) and (30) we get the value of Om. Fig. 13 represents the variation of Om against H. The value 
of parameters are taken as, a = .5, Wm = -01, n = .5, Pcard = .1, A = .01, po = .001, k = 1,1 = .9, Hg = 72. From 
the figures, we see that g, Wc, w, q and Om decrease as H decreases and lies between and 1. The parameter 
s decreases and keeps positive sign as r decreases from positive to negative values. 



C. Exponential Model: 



In this model g{p) is defined as [26], 



gip) = Pm exp 



(38) 
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where pcard is a characteristic constant energy density and n is a dimensionless positive constants. So from (9), 
(14), (16) and (38) we get, 



and 



where, 



Pe = (-l + e^"-)(l + ^)'('+'"'"Vo 

Pc = -(1 + zf^^^-^^^ pQ{Wra + e^^-''(-W™ + nXE^p{l + Wra))) 
■(l + 2)-3(l+--)p,,,/ 



(39) 
(40) 



X 



Exp 



Po 



So the equation of state is given by, 

-(1 + z)^^(i+-™)po(w,„ + e^^-(-«;„, + nXsrcpil + w^))) 



Wc 



(-l + e^«-j')(l + z)3(i+«'")po 



(41) 



and 



powoil + z)3(i+--) - (1 + 2)3(i+--)po(w;m + e^^-(-«;^ + nXs,p(l + u.„0)) 



Po(l + ;2)3(i+'"™) + (-1 + e^«-p)(l + z)3(i+'"™)po 
From (16), (18), (19), (20), (39) and (40) we get the statcfinder parameters and deceleration parameter 



(42) 



and 



r- = - (2 + 9^™ + 9(w^ + n^Xl^pil + w„,f + nXsxpil + w,n){-l + n + (-2 + n)w,„))) (43) 



S = i (ll ~9n + 9Wm - QnWm - 9nXExpi^ + Ifm) 9"-^>n(l + ^rn) 



-Wm + nXExpi'^ + Wjn) 



q= - [1+ 3Wrn - 3nX_Ea;p(l 



p\-L -r uJmJ 



(44) 
(45) 




FIG. 14: The variation of w against H for Wm ~ .01, n = 2, pcard = .1, A = .01, po = .001, k = 1,1 = .9, Ho — 72. 

Fig. 14, 15 and 16 represents the variation of 5, Wc and w against H, z and z respectively. From (9), (10), 
(16), (39) and (40) we plot the graph which shows the variation of the square of velocity of sound against 
H is given in Fig. 17. Fig. 18 represents the variation of s against r and Fig. 19 represents the variation of q 
against H. From (8), (21) and (38) we get the value of Om. Fig. 20 represents the variation of Om against 
H. The values are taken as, Wm = .01, n = 2,pcard = .1,A = .01, po = .001, fc = IJ = .9,Ho = 72. From the 
figures, we see that g,Wc,'w,q and Om decrease as H decreases and lies between and 1. The parameter s 
increases and keeps negative sign as r decreases. 
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FIG. 15: The variation of w against z for w,n = .01, n — 2, pcard = .1, A = .01, po = .001, fc = 1, Z = .9, Ho = 72. 
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FIG. 16: The variation of w against z for Wm = .01, n = 2, pcard = .1, A = .01, po = .001, fc = 1, Z = .9, -ffo = 72. 




FIG. 17: The variation of against H for Wm = .01, n = 2, p^ard = .1, A = .01, po = .001, fc = 1, Z = .9, i/o = 72. 
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FIG. 18: The variation of s against r for Wm = .01, n = 2, pcard = .1, A = .01, po = .001, A: = 1, Z = .9, Hq = 72. 
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FIG. 19: The variation of q against z for Wm = .01, n = 2, pcard = .1, A = .01, po = .001, fc = 1, Z = .9, -ffo = 72. 
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FIG. 20: The variation of Om against H for Wm = .01, n = 2, p^ard = .1, A = .01, pa = .001, fc = 1, Z = .9, i^o = 72. 



D. Modified Exponential Model 

111 this model g{p) is defined as [26], 

9{P) = {Pm + Pcard) exp ( ^P"'^'''^ j 

. \ Pm + Pcard J 



(46) 



where Pcard is a characteristic constant energy density and a and n are dimensionless positive constants. So 
from (9), (14), (16) and (46) we get. 



Pc - -(1 + zf^^+^-^po + eV<i+^)^<^+™"Vo+._.; ((1 + + p,,,,) 



(47) 



and 



= (1 + z)3-™(-(l + zfp^W^ - e^^-^^((l + z)-3"''"pearrf + (1 + z)'po(-U.„ + nXM£)"(l + W,n)))) (48) 

So the equation of state is given by, 

(1 + z)3--(-(l + z)^poWm - €^"^{{1 + z)-3'"" Pcard + (1 + z)^ po{-W,n + uXMsil + W„)))) 



Wc 



-(1 + z)3(i+-™)po + e^M^((l + z)3(i+-^)po + p,,,^) 



(49) 



and 



where 



,Po(l + 2)3(1+-"-) ^ (1 ^ 2)3-^(_(i + zfpoW„, - e^--((l + z)-3--pc„,d + (1 + zfpo{-w.,n + nXuE^X + 



P0(1 + Z)3(l+-..) - (1 + z)3(l+-™)po + e^M«((l + z)3(l+-".)po + Pcard) 



(50) 



MB 



^Pcard 



(l + z)3(l + ^"Opo + p,„,rf 
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FIG. 25: The variation of s against r for a = 2.5, Wm = .01, n = .8, pcard = .1, A = .01, po = .001, k — 1,1 — .9, Ho = 72. 




FIG. 26: The variation of q against for a = 2.5, Wm = .01, n = .8, pcard = .1, A = .01, po = .001, = 1, / = .9, ffo = 72. 
From (16), (18), (19), (20), (47) and (48) we get the statefinder parameters and deceleration parameter, 

r = [2(1 + z)"'''"'"pLrrf + (1 + 2)'""""P0Pcard(4 " 9(-l + nXME)w^{l + W^)) + (1 + zfpl{2 + 9Wra 

+9(w^ + n^Xlsi'^ + w„,f + nXME{l + w„,){-l + n + {-2 + n)w„,)))] / [2(1 + 2)-^"''"pLrd + 2(1 + zfpo] 

(51) 

and 

S = [2(1 + z)-^"''"pLrd + (1 + ^)'"^"'"'P0Pcard(2 + nXME{l + Wr,^){2 + 9Wr,^) - W™(11 + 9w„)) 

-(1 + z)Vo(9n^^MB(l + w^m)^ + Wm{n + 9w,n) + TlXMsil + W™)(-11 + 9n + 9(-2 + n)w,n))] 

I [(2((1 + zfp^ + (1 + z)-''"'"Pcard)((l + ^)"''""Pcard + (1 + zf Pq{-W„^ + nXMf;(l + ^m))))] (52) 



2(1 + z)~^^'"Pcard + (1 + zfpQ{-\ - 3u>,„ + 3nXMg(l + w„0) 
2(1 + Z)3p0 + 2(1 + z)-3»'"p,ard 



(53) 
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FIG. 27: The variation of Om against H for a = 2.5, ~ .01, n = .8, Pcard = .1, A = .01, po = .001, fc = 1, Z = .9, i/o 
72. 
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Fig. 21, 22 and 23 represents the variation of g, Wc and w against H respectively. From (9), (10), (16), (47) 
and (48) we plot the graph which shows the variation of the square of velocity of sound against H is given 
in Fig. 24, Fig. 25 represents the variation of s against r and Fig. 26 represents the variation of q against H. 
From (8), (21) and (46) we get the value of Om. Fig. 27 represents the variation of Om against H. The values 
are taken as, a = 2.5, Wm = .01, n — .8, Pcard = .1, A = .01, po = .001, fc = 1, / = .9, Hq = 72. From the figures, 
we see that g, Wc, w, q first decrease then increase and Om decreases as H decreases and lies between and 
1. The parameter s increases and keeps positive sign as r decreases. 



V. DISCUSSIONS 

In this work, the we have considered Cardassian Universe in Hofava-Lifshitz gravity. The energy density and 
pressure for Cardassian term have been found. Four types of Cardassian Universe models i.e., polytropic/power 
law, modified polytropic, exponential and modified exponential models have been considered for accelerating 
models. To investigate the natures of statefinder parameters, deceleration parameter, Om diagnostic and EoS 
parameters for all types of Cardassian models in Hofava-Lifshitz gravity, we have drawn all parameters w.r.t. 
Hubble parameters H. 

In polytropic/power law Cardassian model of the universe g^w^q and Om parameters have been drawn in 
figures 1, 2, 5 and 6. We have seen that g,w,q and Om decrease as H decreases. The values are taken as, 
Wm = .01, n — .5, Pcard = .1,A = .01, po ~ .001, fc = 1,1 = .9, ifg = 72. In this model, the EoS parameter Wc 
is constant. The variation of the square of velocity of sound against H is given in fig. 3 and it has been 
observed that lies between and 1. Also fig. 4 represents the variation of s against r and the figure shows 
that s increases and keeps positive sign as r decreases from positive to negative values. 

In modified polytropic Cardassian model of the universe g, Wc, w, q and Om parameters have been drawn in 
figures 7, 8, 9, 12 and 13. We have seen that g,Wc,w,q and Om decrease as H decreases. The values are taken 
as, a — .b,Wm — .01, n = -^jPcard = .1,A = .01, Po ~ .001, fc = 1,/ = .9, iJo ~ 72. The variation of the square 
of velocity of sound against H is given in fig. 10 and it has been observed that Cg lies between and 1. Also 
fig. 11 represents the variation of s against r and the figure shows that s increases and keeps positive sign as r 
decreases from positive to negative values. 

In exponential Cardassian model of the universe g,Wc,w,q and Om parameters have been drawn in figures 
14, 15, 16, 19 and 20. We have seen that g,Wc,w,q and Om decrease as H decreases. The values are taken as, 
Wm = .01, n = 2, Pcard = .1,A = .01, po — .001, fc ^ 1,1 = .9, i^Q = 72. The variation of the square of velocity 
of sound Cg against H is given in fig. 17 and it has been observed that lies between and 1. Also fig. 18 
represents the variation of s against r and the figure shows that s increases and keeps negative sign as r decreases. 

In modified exponential Cardassian model of the universe g,Wc,w,q and Om parameters have been drawn in 
figures 21, 22, 23, 26 and 27. We have seen that g, Wc, w, q first decrease then increase and Om decreases as H de- 
creases. The values arc taken as, a = 2.5, Wm = .01, n = .8, Pcard = .1, A = .01, po = .001, fc = 1, / = .9, Hq = 72. 
The variation of the square of velocity of sound against H is given in fig. 24 and it has been observed that cj. 
lies between and 1. Also fig. 25 represents the variation of s against r and the figure shows that s increases 
and keeps positive sign as r decreases. 
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